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The excessiveness of integration-by-part (IBP) identities is discussed. The Lie-algebraic structure 
of the IBP identities is used to reduce the number of the IBP equations to be considered. It is shown 
that Lorentz-invariance (LI) identities do not bring any information additional to that contained in 
the IBP identities, and therefore, can be discarded. 
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O ■ I. INTRODUCTION 

o : 

Calculation of the multiloop radiative corrections in different physical processes becomes more and more important 
bJQ, nowadays. Mainly this is because of the increasing precision of the modern experiments, both in high-energy physics 
and in spectroscopy. The use of the IBP identities 0, Q is a standard approach to the effective calculation of the loop 
\ integrals. These identities can reduce the problem of calculation of arbitrary integral with a given topology to that 
of calculation of the limited number of simpler integrals of the same topology and its subtopologies. 

However, the application of the IBP identities is hampered by their infinite number. The problem is that it is not 
i— i always clear which identities should be used to reduce a given integral. One standard approach, which has proved 
' to be useful, is considering the identities starting from the simplest ones and creating a database of the rules for the 
reduction [3(. This algorithm is essentially sequential, since, in order to solve the next identity, it is necessary first 
O i' to substitute all integrals which are already in the database. Another approach to the problem is to reduce it to the 
problem of "division with the remainder" with respect to some ideal [j, [a, @ . The main difficulty on this way is to 
derive the Grobner basis of the ideal allowing for the correct determination of the simplest remainder. According 
to Ref. 0], this procedure essentially depends on the ordering chosen. This choice is, in some cases, not a simple 
£Nl 1 problem, and the Grobner basis then can be hardly found. Thus, the first approach (Laporta approach) appears to 
be necessary, at least, in this situation. 

One of the problems which can spoil the effectiveness of the Laporta approach is the excessiveness of the IBP 
identities. When the IBP identities are considered one-by-one, most identities do not give any additional information. 
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Indeed, the number of the identities grows as the volume of the region in 1 N multiplied by the number of identities in 
one point, L(L + E) (L being the number of loops, E, the number of external momenta), while the number of integrals 
involved grows only as the volume of the region in 1 N . Thus, in the asymptotics, only one identity out of L(L + E) 
give new information. Other identities in the Laporta approach are checked to reduce to = 0. Unfortunately, this 
check can be a very time consuming calculation. The determination of the minimal set of the IBP identities is also 
• • , important for the analytical solution of identities (i.e., the derivation of the reduction rules) as it allows one to consider 
. ^ ' a smaller set of the identities. 

Though the algebraic manipulation with the IBP identities has been used for a long time, the observation that IBP 
identities form a closed Lie algebra has not been exploited so far. The main purpose of this paper is to demonstrate 
how the Lie-algebraic properties of the IBP equations can be used for both "division with the remainder" algorithms 
and Laporta-like algorithms. 

II. GENERAL DISCUSSION 

Assume that we are interested in the calculation of the L-loop integral depending on the E external momenta 
Pi, . . . ,Pe- There are N scalar products depending on the loop momenta If. 

Sik = h ■ Ik , 1 < i < L, k < L + E, 

N = L(L + l)/2 + LE (1) 
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where gi,...,L = qL+i,...,L+E = Pi,...,E- 

The loop integral has the form 



J(n) = J(n u n 2 ,. .. ,n N ) = / d v h . . .d v l L j{n) = I * n " ' — (2) 



where the scalar functions D a are linear polynomials with respect to Sy . The functions D a are assumed to be linearly 
independent and to form a complete basis in the sense that any non-zero linear combination of them depends on the 
loop momenta, and any can be expressed in terms of D a . Thus, each integral is associated with a point in W, N . 
Some of the functions D a correspond to the denominators of the propagators, the other correspond to the irreducible 
numerators. E.g., the iif-legged L-loop diagram corresponds to E — K — 1 and the maximal number of denominators 
is M = E + 3L — 2, so that the rest N — M = (L — 1)(L + 2E — 4)/2 functions correspond to irreducible numerators. 
For vacuum diagrams, M = 3(X - 1), and N — M = (L— 2){L - 3)/2. 

The IBP identities are based on the fact that, in the dimensional regularization, the integral of the total derivative 
is zero. They are derived from the identity 

= J d v h . . . d v l L O lk j(n) = J d v h . . . d V l L -^- ■ q k j(n) . (3) 

Performing the differentiation in the right-hand side and expressing the scalar products via D a , we obtain the 
recurrence relation for the function J. 

There is also another class of identities, called Lorentz-invariance (LI) identities due to the fact that the integral 
([2]) is Lorentz scalar [8J . They have the form 

PiPj \ ^2,Pk[v—fi ) J(ni,n 2 , ■ • -,n N ) = (4) 

The differential operator in braces is nothing but the generator of the Lorentz transformation in the linear space 
of scalar functions depending on pk- If we explicitly act by the differential operator on the integrand, we obtain 
LI identity. Though these identities can be convenient in some cases, they can be easily represented as some linear 
combination of the IBP identities (see Appendix) and will not be considered in the following. 

For the reduction procedure to work, it is necessary to define some suitable ordering of the integrals, i.e., the 
ordering in Z N . First, one introduces the notion of sectors in Z . The (Ox, ■ ■ ■ , 9n) sector, where &i = 0, 1, is a set of 
all points (ni, . . . , n^r) in Z N whose coordinates obey the condition 

sign (n a - 1/2) = 20 a -l. (5) 

In particular, the point [Q\, . . . , 9n) belongs to the (6\, . . . , 9n) sector, and can be referred to as the corner point 
of the sector. Owing to this definition, the integrals of the same sector have the same number of denominators. It is 
natural to consider the integrals with less denominators to be simpler. When the number of denominators coincides, 
we will consider the integrals with smaller total power of the numerators and denominators to be simpler. Then goes 
the number of the numerators and the last is the lexicographical ordering. Thus, two points n = [n\ 1 n2 1 ■ ■ ■) and 
n' = (n[, n' 2 , . . .), are said to be ordered as n -< n' iff there exists «o, — 2 < io ^ N, such that n, < n^and for any 
i, —2 ^ i < io holds ni = n\. Here n_2, Ji— i, and Uq are determined as 

N N N N 

n_ 2 = ^ 6 (n - 1/2) = 0*, n^ = J2\n a \, n = ^ 6 (-n a + 1/2) . (6) 

a—l a—1 a—1 a — 1 

The integral J (n) is considered to be simpler than J (n 1 ) if n ~< n' . According to this ordering, the integral 
J(6>i, . . . , 6jy) is the simplest integral of (0i, . . . , 8n) sector. 

III. OPERATOR REPRESENTATION 

Let us introduce, similar to Ref. the operators A a and B a acting on functions in Z N as follows 

(A a f) (ni, . . . , n N ) = n a f (m, . . . , n a + 1, . . . , n N ) , 

(B a f) (m, . . . , n N ) = f (m, . . . , n a - 1, . . . , n N ) . (7) 
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Note that these operators act on function, but not on its arguments, and should not be confused with the conventional 
index shifting operators. Using these operators, we can express the IBP identities as constraints on the function 
J having the form 

-P. J = 0, 

P = a al} A a B,3 + b a A a +c, 

where a a/3 , b a , c are some coefficients. We will denote the operator, corresponding to the Oik as Pit'- 

- (P lk J) (n) = J d v h . . . d v l L O lk j (n) . (8) 

Note that the operators A a , B a form Weyl algebra, 

[A a ,B p ]=S a i3. (9) 
Let C be the left ideal generated by operators Pj&, i.e. a set, consisting of all operators, which can be represented 



as 



J2°ikPik, (10) 

where Cik are some polynomials of Ax,.. ■An,B\ 1 . . . ,Bn- This ideal has a simple meaning: for any L G C the 
relation 

(LJ){n 1 ,...,n N ) = (11) 

is a linear combination of some IBP identities. In fact, any linear combination of the IBP identities can be represented 
in a more specific form 

(LJ)(1,...,1) = 0, (12) 

since shifting of the indices can be done by acting from the left with some powers of A a or B a . At first glance, the 
problem of reduction is equivalent to that of division with the remainder by the ideal £, which is effectively solved 
by the construction of the Grobner basis. However, there is an additional obstacle. Note that for any function / of 
N integer variables the following relation holds 

(B a A a f) (1, . . . , 1) = . (no summation) (13) 

Indeed, 

(B a A a f)(l,...,l) = (A a f)(l,...,0,...,l) =0 x/(l,.. .,!,.. .,l) -0, (14) 

where the overscript a denotes the position of the index. Let 1Z be the right ideal generated by the elements 
(BiAi) , . . . , (BnAn). By definition, it consists of all operators of the form 



R = Y,B a A a C cn (15) 



where C a are some polynomials of A\, . . . An, Bi, . . . , Bn. It follows from Eq. ([13]) that 

(i?/)(l,...,l) = 0. (16) 

Thus, for the reduction procedure to work, we have to have an algorithm of division with the remainder by the 
direct sum of the left ideal C and the right ideal 1Z. That means that we have to invent the algorithm allowing the 
decomposition 

p = L + R + r, (17) 

where L G C, R G 1Z, and r is the simplest possible with respect to the ordering chosen. Even though the problem is 
clearly formulated, such algorithm appears to be unknown so far. 
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IV. LIE-ALGEBRAIC STRUCTURE OF THE IBP IDENTITIES 

The operators 

ik = J- • ft (18) 
oh 

form a closed algebra with the commutation relations 

[O ik ,O j i] = 5 a O jk -8 jk O il . (19) 

We can easily check that the operators Pi k obey the same commutation relations as Oik ■ This algebra is nothing but 
the algebra of the group of linear changes of variables 

h -> M lk q k . (20) 

The operator Oi k corresponds to the infinitesimal transformation U — > = li + eq k in the sense that 

d 



f(s' lm )d D l' 1 ...d v l' L =if(s l 



dh qkf{Sh * 



d u l l ...d u l L + 0{e 2 ). (21) 



Note that the so-called symmetry relations are also the consequences of the invariance of the integrand under the 
action of some elements of this group. 

The scaleless integral can be defined as the one which gains additional non-unity factor under some transformation 
([2H)l . This definition corresponds to the conventional notion of scaleless integrals. E.g., owing to this definition, the 
integral 

I= f *h*h ( 22) 

J (h-pfllih-p-h) 2 

is scalelless as, under the transformation 

li — > ali + (1 — a ) P, 

h -» al 2 , (23) 

it transforms as 

/ -» a 2V - 6 I. (24) 

In dimensional regularization, the scaleless integrals are zero, as well as the integrals which differ from the scaleless 
ones by additional polynomial factor in the numerator. Thus, once the integral in the corner point of the sector is 
scaleless, the whole sector is zero. Let us prove a simple criterion of zero sectors. 

Criterion 1. If the solution of all IBP relations in the corner point of the sector (Qi,..., On) results in the identity 

J(6 U ...,9 N ) = 0, (25) 
then this sector is zero, i.e., all integrals of this sector are zero. 

Indeed, by the condition, j(0±, . . . ,6n) can be represented as the action of some linear combination of Oi k on 
j(9\, . . . ,8n). Since these operators are generators of the transformation (|20|) . we can conclude that J{0\, . . . , 0jv) is 
scaleless and thus, the whole sector (9%, . . . , 9n) is zero. This criterion gives a simple and convenient way to determine 
zero sectors. 



V. EXCESSIVENESS OF THE IBP IDENTITIES SET. 



In this Section we describe some consequences of the algebraic structure of the IBP identities. 
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Proposition 1. Let L ^ 2. Then of all L(L + E) IBP identities we can consider only identities, generated by the 
operators: 

d 

— — ■ h+i, i = l,...,L, Il+i = Ii 
ok 

t^-k (26) 

i— 1 

Indeed, this set of operators form the multiplicative basis of the Lie-algebra ([!§]) , i.e., the rest of the operators 
can be obtained from the commutators of the chosen ones. The total number of the operators in the set (|26|) is 
L + E + 1, which is smaller than the original L{L + E) for L 2. This simple fact can be used for the construction 
of the reduction rules and also for the Laporta algorithm. Nevertheless, such system of the IBP identities is still 
overdetermined. In the asymptotics only 1/(L + E + 1) part of the identities gives new information. 

Now we prove a more refined criterion for the identities which can be thrown away without loss of information. Let 
S = {Pi, Pi, . . . , Pk} be some set of the IBP operators, and P be some IBP operator with the following property: its 
commutator with any Pk G S is a linear combination of Pi G S. Then we have the following 

Criterion 2. If for some point n 6 1, N the integral J(n) can be expressed via simpler integrals with the help of the 
identities obtained from the operators in S, then the identity (PJ)(n) — can be represented as a linear combination 
of the identities obtained from the operators in S and the identities of the form {PJ){n') — with n' -< n. 

Proof. To prove it, we note that, since the integral J(n) can be expressed via simpler integrals with the help of the 
identities obtained from the operators in S, there exist such polynomials Ci that for any function / 

J2(C l P l f)(l) = f(n) + o(n), (27) 
i 

o (n) denotes here the linear combination of f(n') with n 1 -< n. Now we substitute / = PJ: 

J2(C l P l PJ)(l) = PJ(n)+Po(n), (28) 
i 

and use the commutation relation 

P l P = PP l +Y,ci m P m (29) 

m 

where ci m are some constants. We obtain 

PJ(n) = Po (n) + 53(C,Pfl J)(l) + J2(C lClm P m J){l) (30) 

l Ira 

Here Po (n) denotes the linear combination of PJ(n') with n! -< n, and the last two terms is a linear combination of 
of the identities obtained from the operators in S. □ 

Let us consider the sequence of the operators 

{Pi,... P N } = {Pl,L+E, ■ • • , Pl.l,P2.L+E, • ■ • , P2,2, P L .L+E, P L , L } (31) 

Note that the number of the operators in this sequence equals to N = L(L + l)/2 + LE, the total number of the 
denominators and numerators in basis. This sequence has the following property: for any Vi the criterion [3] applies 
with S = {Vi, . ■ .Vi-i}. Suppose we have solved the identities, generated by operators {Vi, . . .Vi-i}, then the 
identities, generated by Vi should be solved only in points for which there are no reduction rules yet. On each step 
of this procedure the "dimension" of the set of such points is decreased by one, thus, when we have considered all 
operators in the sequence, we have only finite number of the integrals, which are not yet reduced. The rest L(L — l)/2 
IBP identities can be used for the reduction of integrals in these points. Note, that the mean number of the identities, 
considered in each point is 1, as it should be. The choice of the sequence (|31l) is, of course, not unique. 

Now we derive the criterion which can be used for the algorithm combining the "division with the remainder " and 
Laporta method. 
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Criterion 3. Let in some sector the identities generated by some operator P — CikPik have the form 

(PJ){n) = J(n) +o(n) (32) 

and for any n in the sector holds 

n -< h. (33) 

Then the identity generated by another operator P' — ^2c' ik Pik ^ P in point h, corresponding to the integral 
expressed by P , can be represented as a linear combination of the identities of the form PJ and the identities in 
simpler points. 

Proof. To prove it, let us consider the identity 

{PP'J) (n) = (P'J) (n) + (P'o) (n) = P'PJ (n) + ([P, P'] J) (n) . (34) 

Thus 

(P'J) (n) = {P'PJ) (n) - (P'o) (n) + ([P, P'} J) (n) . (35) 

The first term in the left-hand side is some linear combination of the IBP identities generated by P, and the last two 
terms contain only identities in the points n' -< n. Thus, the identity P'J (n) is dependent on the identities of the 
form PJ and on the identities in simpler points. □ 

The basic idea of application of this criterion is the following. Consider the IBP identities in general point of the 
sector. We can either use the set of identities, generated by Pik or by other L(L + E) independent linear combinations 
of Pik- In fact, it is natural to pass to such linear combinations, in which all most complex integrals are different 
("solve" the identities in general point). Among these identities, select one, having the form (|32| with n -< n and 
fi corresponding to the integral not yet expressed. Determine the points of the sector corresponding to the integrals 
for which this identity does not work. Consider the other identities only in these points. Repeate the same steps, 
starting from the selection of the identity. After some iterations, we might be unable to find the identity matching the 
conditions. At this stage, we have to solve the rest of the identities only in the points, corresponding to the integrals, 
which are not yet expressed. E.g., we can use the Laporta method. The advantage of this combined approach is that 
the "dimension" of the set of points in which we use the Laporta method is usually essentially less than N. 



VI. CONCLUSION 



In the present paper, we have considered the dependencies between the IBP identities. Account of these dependencies 
dramatically decreases the number of the identities to be considered. They come from the fact that the corresponding 
operators Pik form a closed Lie algebra of the group of linear change of variables. Using this interpretation of the 
IBP identities, we have proved the simple criterion of the zero sectors. The two criteria of the excessiveness of an 
identity were proven. Probably, using the Criterion [2l it is possible to prove the finiteness of the number of master 
integrals in general case. Indeed, selecting the sequence of the operators as it was described, we decrease on each step 
the "dimension" of the set of the unexpressed integrals by one. Since the number of the operators in this sequence 
equals N, we are left with the set of "dimension" zero, i.e., consisting of finite number of points. The Criterion [3] 
can be used for the combined algorithms in which the Laporta reduction is performed on some subset of the points 
of the given sector with the dimension less than N. It was also shown that the Lorentz-invariance identities can be 
completely discarded. 

This work was supported by RFBR Grant No. 07-02-00953. I also thank for warm hospitality the Max-Planck 
Institute for Quantum Optics, Garching, where a part of this work was done. 

APPENDIX A: EXPRESSING LI IDENTITIES VIA IBP IDENTITIES. 

In this Appendix we will show that the LI identities always can be represented as a linear combination of the IBP 
identities. 

Let us note that the integrand j in Eq. {2J is a scalar function of qi. Thus, the operator 

L+E d L d E d 

H^TT] =E^7^] ( A1 ) 

k=l dq£ fe=1 dl^ fe=1 Op^ 
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when acting on the integrand j, annihilates it identically. Indeed, this operator is nothing, but the generator of 
Lorentz transformations in the linear space of functions of Thus, the operator in Eq. (|H), when acting on the 
integrand, can be represented as follows: 



nv 3 



^3 L-\-E q 



k=l 



d Pk 



fe=i 

L 



k=l 



dq 



k=l 



01 



E 

L 

E 

fc=i 



(Pi ■ h)Pj 







Oh 



(pj -h) Pi 







dl k 



■ Pj {Pi ■ h) - — ■ Pi (Pj ■ h) 



Oh 



Taking into account that the scalar products (pij - Ik) in the last line can be expressed via D a , we conclude, that 
any LI identity can be represented as a linear combination of the IBP identities. 
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